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Quenched QCD with domain-wall fermions on coarse lattices * 
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We investigate the existence of chiral zero modes at a -1 ~ 1 GeV in quenched domain-wall QCD. Simulations 
are carried out for the plaquette and an RG-improved gauge actions on a 12 3 x 24 x N s lattice with N s = 10 - 50. 
We find that the pion mass in the chiral limit remains non-vanishing as iV s — » oo for both gauge actions. Possible 
origins of this non-vanishing pion mass are discussed. 



1. INTRODUCTION 

The domain-wall fermion formulation of QCD 
(DWQCD) jl]j|] is expected to realize exact chiral 
symmetry on the lattice without species doubling 
at finite lattice spacing. This represents an ap- 
pealing possibility, particularly for investigations 
of problems such as weak matrix elements sensi- 
tive to chiral symmetry || . Therefore, a number 
of studies have been made Q . 

Simulations in DWQCD, however, requires 
considerable computing power. Even if the size 
of the extra dimension N s can be taken as small, 
i.e. N s = O(10), lattice spacings much finer than 
a -1 w 2 GeV will be difficult to simulate even 
in quenched QCD. Hence simulations on coarse 
lattices down to a -1 k 1 GeV will be needed for 
reliable continuum extrapolations. 

A first step in DWQCD at such a strong cou- 
pling will be to ensure the existence of chiral zero 
modes. Here we report our preliminary results on 
this problem. Our study is made for the plaquette 
action, and also for an RG-improved action || to 
examine the effects of reduced discretization er- 
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2. PARAMETERS 

We employ the fermion action identical to that 
in Ref. g], with the domain wall height M and 
bare quark mass to/. The gauge coupling is cho- 
sen to be [3 — 5.65 for the plaquette action and 
= 2.2 for the RG action; these values corre- 
spond to the scale a -1 « 1 GeV determined from 
the string tension. Simulations are made on an 
12 3 x 24 x N s lattice with N s = 10, 20, 30 and 50. 

In the free fermion case the chiral zero mode ex- 
ists over the range < M < 2. Since this range 
will be shifted in the interacting theory, we em- 
ploy M = 1.3, 1.7, 2.1 and 2.5. For each of these 
values of M, we take mj = 0.1, 0.05 and 0.03 and 
calculate the pion mass for both degenerate and 
non-degenerate quark and antiquark pairs. For 
each parameter point we have typically 20 config- 
urations. 



3. RESULTS 

In Fig. [l] the pion mass squared to^. is plotted as 
a function of the averaged bare quark mass mj" 
at M — 1.7 in the case of the plaquette gauge 
action. Since the linearity of in mj" is well 
satisfied, we adopt a linear chiral extrapolation in 
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Figure 1. Pion mass squared as a function of m^" 
at M = 1.7 for the plaquette action. Solid circles 
show the N s —> oo limit of open symbols for each 
rrif. Lines show the linear fits. 
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As it can be seen in the Fig. [j], we find a 



non-zero value for m'i at 
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for each N s 



This non-zero value, however, is expected to 
vanish exponentially as N s — » oo. To see whether 
this occurs, we plot mi at mf — by thick 
square symbols in Fig. ||(a) as a function of N s . 
Fitting these points by the form ae~^ Ns yields 
the dotted line with x 2 /dof = 21.4. An alterna- 
tive fit, allowing a constant, c + ae~^ Ns gives the 
solid line with y?/dof = 1.7 and c = 0.0532(75). 
Clearly the latter fit better reproduces the behav- 
ior of our data. A similar phenomenon has been 
previously reported also at (3 — 5.7 [||. 

In order to confirm the existence of a non-zero 
c, we attempt to interchange the order of the lim- 
its ui°j v — > and N s — > oo. As shown by solid 
lines going through open circles in Fig. |^(a), we 
first make a fit of form rn%(mj", N s ) = c'(mj") + 
ae - ^' to make an N s — * oo extrapolation for 
each value of mj". The results, shown by solid cir- 
cles in Fig. [l], exhibits a linear behavior in in J"- A 
linear chiral extrapolation rn^{my,N s = oo) = 
d{mf) = d + jmf then yields d = 0.0531(70), 
which agrees well with the value of c previously 
obtained. The commutativity of the two limits 
is summarized in Fig. ||(a), where the two values 
obtained for rrv^(mj v = 0, N s = oo) are marked 
by "X" and "Y". 

The above analyses strongly support the con- 




Figure 2. Pion mass squared as a function of N s 
at M — 1.7 for the plaquette (a) and RG (b) 



actions. Two limits 
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and N„ 



and 



their combinations in two possible orders denoted 
by "X" and "Y" are shown. 



elusion that at m^" = does not vanish in 
the limit N s = oo at a" 1 ~ 1 GeV for the 4- 
dimcnsional lattice size of 12 3 x 24. We find this 
conclusion to apply as well to the RG-improved 
gauge action, as shown in Fig. ||(b). 

In Fig. H the residual m^{jnf — 0, N s = oo) is 
plotted as a function of M for the plaquette and 
RG-improved gauge actions. We observe that a 
non-vanishing residual remains under the varia- 
tion of M for both actions. An improvement may 
be seen for the RG action, however, in that the 
magnitude of the residual is reduced by about a 
factor of two compared to that for the plaquette 
action. 

Finally the decay rate £ extracted from the fits 
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Figure 3. Summary of the squared pion mass in 
the double limit m,f — > and N s — > oo. 

is shown in Fig.^ as a function of m™ at M = 
1.7. For both actions, the chiral extrapolation of £ 
from mj ^ (open circles) smoothly agrees with 
the one directly obtained at m™ = assuming 
c / (filled circle), but disagrees with the value 
from the c = fit. Furthermore the values of £ 
varies little with mf . This is consistent with the 
expectation that the decay rate is governed by 
the transfer matrix in the extra dimension, which 
is independent of mj. 

4. DISCUSSION 

Our study of domain-wall QCD with the pla- 
quette gauge action has shown that the pion mass 
in the chiral limit remains non-zero for an in- 
finitely sized extra dimension at a strong gauge 
coupling corresponding to o" 1 ~ 1 GeV on a 
12 3 x 24 lattice. We have also found that this con- 
clusion remains unchanged for the RG-improved 
gauge action. 

One possibility for the origin of a non-zero pion 
mass is finite spatial size effects, which is be- 
ing checked by increasing the spatial size from 
12 to 16. Another possibility is that it is an 
artifact of the linear chiral extrapolation, which 
does not take into account the possible presence 
of quenched chiral logarithms. To explore this 
possibility, the chiral breaking term in the axial 
Ward-Takahashi identity g, which we expect to 
be free from the quenched singularities, is cur- 
rently being investigated. 

Finally it is possible that no chiral zero modes 
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Figure 4. Decay rate £ from the curves of Figs. |2j 
for the plaquette (a) and RG (b). 

exist at a -1 ~ 1 GeV, and therefore the domain- 
wall formalism fails to realize chiral symmetry in 
the region of strong coupling. A plausible expla- 
nation for this failure might be that the range of 
M for zero modes to exist disappears for stronger 
couplings, where there is no gap of eigenvalues for 
the 4-dimensional Wilson operator Qj. 
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